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ABSTRACT 

Diffraction limited resolution adaptive optics (AO) correction in visible wavelengths requires a high performance 
control. In this paper we investigate infinite impulse response filters that optimize the wavefront correction: 
we tested these algorithms through full numerical simulations of a single-conjugate AO system comprising an 
adaptive secondary mirror with 1127 actuators and a pyramid wavefront sensor (WFS). The actual practicability 
of the algorithms depends on both robustness and knowledge of the real system: errors in the system model may 
even worsen the performance. In particular we checked the robustness of the algorithms in different conditions, 
proving that the proposed method can reject both disturbance and calibration errors. 

Keywords: Adaptive optics, Pyramid wavefront sensor, telescope vibrations, adaptive secondary mirror, modal 
control optimization 

1. INTRODUCTION 

Modern 8m telescopes are able to provide up to 90% Strehl Ratio (SR) on high flux regime in H band (Esposito 
et a?P}, and reach a maximum resolution of about 30mas (FLAO@LBT at the J band or the Hubble Space 
Telescope at the U band or using lucky imaging technique on large telescope at R band Law et alW). Since 
the maximum achievable resolution is fixed by diffraction retrieving Point Spread Function (PSF), which Full 
Width at Half Maximum (FWHM) is X/D, where A is the imaged wavelength and D the telescope diameter. 
8-m class telescope could improve the maximum achievable resolution making images at shorter wavelenghts, 
such as V band, which could provide a 12mas resolution. However behind these numbers is hidden a limitation: 
an Adaptive Optics (AO) system for visible wavelengths has tighter requirements compared to an near infra-red 
one. In fact shorter wavelengths correspond to smaller Fried's parameters, ro, hence more degrees of correction 
and low measurement noise values are required to obtain a good correction. 

Optimized control techniques such as a Modal Control Optimization^ observer-based and model-based^^' 
approaches can increase the performance of an AO system in comparison to an integrator controller. These 
methods need to work in realistic conditions, i.e. without an accurate knowledge of the real system and with 
complex disturbances like calibration errors and telescope vibrations. In particular, Modal Control Optimization 
has only one degree of freedom per mode, that is the integrator gain, so it cannot efficiently reject disturbance 
condition such as telescope vibrations, while the other methods rely on the accuracy of the system model. 
For these reasons, starting from the work of Dessenne et oZ.p' we have chosen to design a method that is data 
driven, so it relies as little as possible on the system model, and can have many degrees of freedoms as to adapt 
itself to the disturbances. 

The remainder of this paper is as follows. Sec{2] introduces the Infinite Impulse Response (IIR) filter based 
control along with other controllers that we will use in the performance evaluation as comparison. Sec|3] briefly 
describes the AO system considered in our simulations. Finally in Sec|4] the controllers are analyzed through 
numerical simulations from both performance and robustness viewpoint. 

2. AO CONTROL SYSTEMS 

The AO system control loop considered in this work is represented in Fig. [T] The controller gets the modal 
measurements as input and gives the modal commands as output. The most used controller in AO systems is 
the integrator, so we have chosen it to make a comparison that will be presented in the next sections. 
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Figure 1. Adaptive Optics Control loop. 



2.1 IIR filter based control 



To determine the optimal control expressed as an IIR filter we realized an algorithm that minimized the residual 
variance for each mode i as a function of the IIR filter parameters: 



[ai,6i] = min Jj , 

[a, A] 

where a, and bi are two vectors of the IIR filter denominator and numerator coefficients: 
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where N a and Nb are the order of respectively the denominator and numerator. When the modes are orthonor- 
mal, minimizing the sum of the residual modal variances is equivalent to minimizing the residual phase variance 
in the direction of the Natural Guide Star (NGS), which is the goal of the control system studied in this paper. 
The cost, that is the measured residual modal variance of the i-th mode Ji, is: 



Ji = J2 'K'" 



M = £ \\W(z)\\Hf(u). 



(3) 



where: 



W(z) 



(4) 



l + C(z)WFS(z)DM(z)' 

to is the discrete frequency, n is the number of points of the discrete signals, T is the sampling period, $™- eas (w) 
is the Power Spectral Density (PSD) of the closed loop i-th modal coefficient determined from the WFS mea- 
surements, and is the PSD of the open loop i-th modal coefficient. The open loop measurement o, would 
be the input signal of the AO system without the control feedback, and we reconstruct it as the sum of the 
closed loop measurement s, and command q. At iteration k, taking into account the delay d, it can be written 
as: 

Oi(fc) = Si (k) + Ci(k - d) . (5) 
It can be shown (see Dessenne et atf^) that: 



Si (k) = W(z) 0i (k) 



(0) 



so that the relationship shown in Eqj3] is demonstrated. 

This algorithm must verify the stability of the loop during the minimization of the cost Ji, and it must discard 
the controller that do not stabilize the loop. In fact, an unstable controller could produce the minimum cost, 
but it could not be implemented. 



To take into account the modelling errors, like pupil shifts or approximations of DM and WFS models, a stability 
constrain is not enough. Therefore the controller must guarantee some robustness, so small variations of the 
model parameters would not drive the controller to instability. Thus we impose max \\fi (Xi( z )) l| 2 < Vt where 
/i(-) gives the roots of the polynomial, < r\ < 1, and Xi( z ) = 1 + Ci(z)P(z) is the characteristic polynomial, 
where P(z) is the AO system model: 

P(z) = WFS(z)DM(z) , (7) 

we have simply assumed that WFS(z) = and DM(z) = z _1 , in order to model a simple AO loop with two 
frames delay. 

Note that -q gives a measurement of the stability margin since if r\ is equal to 1 only stability is guaranteed, 
while smaller values of r\ correspond to a greater distance from instability, i.e. a greater robustness. 
In this work we choose r\ = 0.8, which corresponds to a gain margin of 0.2. Hence the algorithm can be 
summarized as: 



Acquire with a closed loop the mode measurement Sj and commands Cj. 
Determine &f (/). 

• Choose N a and Nj, and the starting IIR filter parameters aj(0), bi(0). 

• Search the combination of IIR filter parameters a.;, bi which minimizes Ji and satisfies max (xi(z)) || 2 < 
V- 

In this work, the minimization of Ji is based on the downhill simplex method of Nelder & Meadf^l and we opted 
for N a — Nb — 3 in case of turbulence and 2 coefficients more for each vibration. So, in the next simulations, 
Tip and Tilt have N a = Nf, = 7 while the other modes have N a = Nb = 3, because, as presented in Secj3j Tip 
and Tilt are affected by two vibrations. 

2.2 Modal Gain Integrator 

The Transfer Function of the modal integrator is: 

C(z) = (8) 
1 — z 1 

where gi is the integrator gain of the z-th mode, and z is the variable of the Z-transform. 

In this work, we refer as integrator the particular case in which all the gains are equal, and as Optimized Modal 
Gain Integrator 3 (OMGI) the particular case in which each modal gain is optimized. For both cases we optimize 
the gain: in the integrator case we optimized the global gain by a trial and error procedure i.e. we ran many 
simulations with different gains and then we chose the one that produced the best results in terms of Strehl 
Ratio (SR); instead, in case of OMGI, we determine the gain following the same algorithm presented in the 
previous section, but instead of minimizing the cost Ji as a function of the IIR filter parameters, we minimized 
it as function of the integrator gain gf. 

(ji = mm Ji , (9) 

9i 

since the integrator controller is an IIR filter with N a — 2 A5 = 1 and a single variable parameter, the gain g. 
An example of IIR and OMGI Rejection Transfer Functions (RTF) W is shown in Fig(2j Note that the OMGI 
give better rejection at the lowest frequencies, but the bandwidth rapidly decreases with the mode order. The 
integrator controller bandwidth is in direct proportion to the peak at high frequencies, and this peak increase 
the residual phase due to the measurement noise. Hence, the high order modes, that are more affected by noise, 
have a lower bandwidth. Instead the IIR filter gives a bandwidth greater than 20 Hz for all the modes with low 
peaks at high frequencies. 




Figure 2. RTFs for various modes with the IIR filter based control with N a = N b = 3 (left), and the OMGI (right). This 
RTFs are determined for the case magnitude 12.5 without vibrations nor pupil shifts. All the other parameters may be 
found in TabE 

3. AO SYSTEM AND SIMULATOR 

The AO system analyzed in this work is based on a pyramid wavefront sensor (WFS) with tilt modulation 
(Ragazzoni^) and on an adaptive secondary mirror (ASM) (Salman et a|Pj). We chose to work with an 8m 
class telescope, an ASM with 1127 actuators, and a WFS with 40 x 40 sub-apertures. As detector of the 
WFS we considered an ANDOR iXon X3 897 camera^ with Read Out Noise (RON), Excess Noise, and Clock 
Induced Charge (CIC) (Daigle et all 1 - 7 -) values shown in Tab[l] All simulations reported hereafter rely on the 
End-to-End numerical simulator used to perform the performance analysis and optimization of the First Light 
Adaptive Optics (FLAO) system of the LBT telescoped 

The atmospheric turbulence was simulated with a set of phase screens with Von-Karman statistics that were 
displaced in front of the telescope pupil to emulate the time-evolving turbulence according to the Taylor hy- 
pothesis. 

Pupil shifts due to the sensor optical misalignments are a common calibration error. We will evaluate in our 
simulations the robustness of the controller to this kind of errors. The pupil shifts considered in the simulations 
have an amplitude of 0.1 pixel in both x and y directions. 

We add in the simulations, as external disturbance, telescope vibrations. Vibrations may arise from many 
different situations, e.g. telescope orientation, telescope tracking errors, and wind shaking. In particular, since 
vibrations cause displacements of the image they have a major impact on so-called tip/tilt modes. Hence, we 
chose to introduce Tip/Tilt vibrations with central frequencies of 13 Hz and 22 Hz, an RMS of 20 mas each, and 
a damping ratio of 0.01. 

4. SIMULATIONS 

Firstly we will check the performance of the integrator with turbulence and without pupil shifts nor vibrations. 
We chose to run all the simulations in between magnitude 8.5 and 12.5. All the simulation parameters are 
listed in Tab{T] The gain is determined with the method described in Sec |2.2| In Fig(3]and Tab[2]are shown the 
results. As we expect the performance of the system rapidly decreases at the shortest wavelengths. Moreover 
the decrease is greater in the visible wavelengths, in fact, while at K band from magnitude 8.5 to 12.5 less than 
8% of SR is lost, at V band the difference is about 40%. 

In the next simulations we considered only these two magnitudes, 8.5 and 12.5, to make a comparison between 
all the controllers. We will compare the controller in 4 difference observation conditions: the first one is with 
turbulence and without pupil shifts nor vibrations, then second and third case with turbulence and with pupil 
shifts or vibrations, respectively, and the last one with all the disturbances included. 

The controllers parameters are optimized for each condition. The gain found for the integrator and the min- 
imum, maximum and mean gain of the OMGI are shown in Tabj3] Regarding the IIR filter Figj4] shown, as 



Telescope 


Effective diameter (D) 


8m 


Central obstruction 


0.138 D 


Pyramid WFS 


Sensing wavelength (A) 


0.75 /jm 


Number of sub- apertures 


40 x 40 


Tilt modulation radius 


2.0A 


Total average transmission 


0.41 


RON (electrons per pixel) 


0.06e"RMS 


Clock Induced Charge 


0.005 e" /pixel/frame 


Excess Noise 


n/2 


Guide star 


Mag. zero point (Johnson^) 


1.7610"* J/s/mV/im 



ASM 


Modes (Karhunen-Loeve) 1127 


Turbulence 


Outer scale (La) 


25 m 


Mean wind speed 


10.6 m/s 


Telescope Vibrat 


on 


Frequencies 


13, 22 Hz 


Standard deviation 


20 mas 


Damping ratio 


0.01 


Loop paramenters 


Sampling frequency 


1000 Hz 


Delay 


2 frames 



Table 1. Summary of the simulation parameters. 
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Figure 3. Summary of simulation results for the integrator controller in case of turbulence without pupil shifts nor 
vibrations. Integrator gain is determined as described in Sec |2.2| Seeing 0.8", see Table [l] for the other parameters. 



an example, two RTFs for the tip/tilt modes in two different cases: on the left, with only turbulence and, 
on the right, with turbulence, pupil shifts, and vibrations. Note that in the first case we design a filter with 
N a = Nb = 3, while with all the disturbances N a = Nf, = 7. 

All the results of the simulations are summarized in Tab[2]in term of SR for all the considered observing con- 
ditions. Fig(5] shows the SR at different wavelengths at magnitude 8.5, Fig(6] shows the corresponding results 
at magnitude 12.5. 

The integrator and the OMGI show similar results in case of turbulence and pupil shifts, while the vibrations 
affect more the integrator. These results is expected, because the OMGI has more degrees of freedoms, one for 
each mode, with respect to the single degree of the integrator. In fact it can be shown in Tab|3]that the OMGI 
maximum gain value is greater in case of vibrations, so it has a greater bandwidth and can better reject the 
Tip/Tilt disturbances. Note that in case of pupil shifts the gains are lower to guarantee greater robustness. 
At the faintest magnitude the differences between the controller are greater, because, as we expected, the IIR 



filters reject better the noise without decreasing the bandwidth (see Sec 2.2 and Fig|2J) . 

The IIR filter based control produces similar performance in all the tested conditions: it can efficiently reject 
both calibration errors and disturbances. In the worst case - turbulence, pupil shifts and vibrations - the IIR 
filter at V band gives the same SR of the OMGI and about two times the one of the integrator both in the 
condition without pupil shifts nor vibrations. 
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94.5 


96.8 


9.5 


50.1 


70.1 


93.7 


96.4 


10.5 


42.4 


64.3 


92.2 


95.5 


11.5 


30.7 


54.4 


89.5 


93.9 


12.5 


14.9 


37.1 


83.3 


90.1 



Table 2. Summary of simulation results for the integrator controller. Integrator gain is determined as described in Sec |2.2| 
Seeing 0.8", see Table [l] for the other parameters. 



Controller 


integrator 


OMGI 


magnitude 


g 


g min 


g max 


g mean 


Turbulence 


8.5 


0.6 


0.23 


0.88 


0.50 


12.5 


0.3 


0.01 


0.85 


0.20 


Turbulence - 


- Pupil shifts 


8.5 


0.6 


0.08 


0.92 


0.33 


12.5 


0.3 


0.03 


0.85 


0.23 


Turbulence + Vibrations 


8.5 


0.8 


0.23 


0.92 


0.50 


12.5 


0.4 


0.01 


0.93 


0.22 


Turbulence + Pupil shifts + Vibration 


8.5 


0.7 


0.08 


0.92 


0.33 


12.5 


0.4 


0.01 


0.92 


0.23 



Table 3. Summary of simulation optimal gain. Comparison between Integrator and OMGI. Seeing 0.8". 



5. PROCESSING POWER ESTIMATE 

In this section we compare the processing power required for the IIR filter based control described above, 
compared to the integrator case. Referencing the scheme in Figjl] we can split the required processing power 
into three major steps: 

1. Reconstruction matrix 

2. Controller 

3. Modes to command projection matrix 

Steps [I] and [3] are common to the different control strategies, while step [2] is where the techniques differ. 
For the estimate, we suppose the following system parameters (partially following Tab[I]): 

• Number of controlled modes: 1127 (including tip-tilt) 

• Number of DM actuators: 1127 

• Number of slopes measured by the WFS: 2480 

• Loop frequency: 1000 Hz 

• System delay: 2 frames 



Transfer function: amplitude Transfer function: amplitude 
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Figure 4. Rejection Transfer Function (RTF) for the Tip mode with an IIR filter based control with N a = Nb = 3 (left), 
and with N a = Nb = 7 (right). Note that the order 7 is used to reject both vibrations at 13 and 22 Hz and atmospheric 
turbulence. These RTFs are determined for the case magnitude 8.5. All the other parameters may be found in Tab[T] 
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Figure 5. Summary of simulation results for turbulence, pupil shifts and vibrations at magnitude 8.5. Comparison 
between Integrator, IIR filter based control, and autogain. Seeing 0.8". 

Step [I] is a straightforward vector/matrix multiplication. The required processing power is iVmodes x -^slopes 
multiply- accumulate (MAC, a basic operation employed by most Digital Signal Processors - DSPs where two 
operands are multiplied and summed to an accumulator). In our example it is about 2.8 Mega-MAC (MMAC). 
Step 3 is also a vector/matrix multiplication, this time of dimension N mo des x -^commands- In °u r example 
this amounts to about 1.3 MMAC. Assuming a total delay of 2 frames at 1 KHz, and further assuming that 
this computation cannot take more than l/10th of the available time in order to be considered negligible, 




Figure 6. Summary of simulation results for turbulence, pupil shifts and vibrations at magnitude 12.5. Comparison 
between Integrator, IIR filter based control, and autogain. Seeing 0.8". 
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12.5 
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8.5 
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Table 4. Summary of simulation results. Comparison between Integrator, IIR filter based control, and OMGI. Seeing 
0.8". 



the processing power requirements for steps [l] and [3] combined is 20 GMAC/sec. This is comparable with 
existing systems, like the LBT one^l where the total processing power, distributed over 168 DSPs, is about 
25 GMAC/sec (only a fraction of this processing power is used on the LBT system, since the reconstrucion 
matrices have smaller dimensions and the DSPs also have more tasks to attend to) . Both steps [I] and [3] can 
be efficiently parallelized over such a large number of DSPs, since each row of the matrix multiplication can be 



processed independently of the others. 

For step [2j the integrator case needs to sum the measurement vector from the previous frame to the current 
one. Additionally, a gain value must be multiplied to each measurement value. This requires about 1 Kilo-MAC 
(KMAC), which is completely negligible (about l/4000th) when compared to the previous numbers. The only 
effect (and a rather large one) is to prevent full pipelining of the two matrix multiplications. 
For the IIR filters, we take into account a filter whith order 7 for the tip-tilt modes, and 3 for the other modes (as 
in Sec |2.1| . With N a — 7 and Nf, = 7, this requires only a few dozens MACs for the tip-tilt IIR filters, and about 
6 KMAC for the other modes. While this number is 6 times bigger than the integrator, it is again negligible 
in the context of the processing power required by steps [I] and [3] Furthermore, no additional requirements are 
posed on the other processing steps, except for the full separation of the two matrix multiplications, which was 
already required for the integrator with modal gain. 

We can therefore conclude that the IIR filtering strategy outlined in this article can readily be implemented on 
existing or future AO systems with little to no penalty in terms of loop delay. 

6. CONCLUSIONS 

In this paper we have shown that IIR filter based control is a control solution that can deliver good performance 
even in difficult conditions. It is a possible choice for AO systems requiring high performance and robustness 
in case of calibration errors and telescope vibrations. 

Moreover, thanks to the data driven approach, the presented method need only a basic knowledge of the system, 
and the computational burden is focused in the parameters optimization that can be made off-line. In fact, the 
computational power requested to the Real Time Computer of the AO system will be negligible in comparison 
to the integrator as we proved in Sec[5] 

As further work, more investigation on the robustness, on the coefficient r\ and on the impact of IIR filter order 
N a , Nb on the performance must be done. 
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